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Abstract 

The phonon dispersion and the electron-phonon interaction for the /3-Po 
and the bcc high pressure phases of tellurium are computed with density- 
functional perturbation theory. Our calculations reproduce and explain the 
experimentally observed pressure dependence of the superconducting critical 
temperature (T c ) and confirm the connection between the jump in T c and 
the structural phase transition. The phonon contribution to the free energy is 
shown to be responsible for the difference in the structural transition pressure 
observed in low and room temperature experiments. 
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High pressure /3-Po and bcc phases of tellurium (Te) are metallic and superconducting 
at low temperatures |1],|2] . The rhombohedral /3-Po structure has one atom per unit cell. 
The bcc structure is a special case of the /3-Po structure with the angle between the lattice 
vectors a = 109.47°. The transition from the (3-Po to the higher pressure bcc phase is 
observed at 27 GPa in X-ray diffraction experiments || at room temperature. At the 
transition the unit cell volume decreases by 2% and the angle a increases by 4°. A jump 
in the superconducting critical temperature (T c ) from 2.5°K to 7.4°K is observed at 32- 
35 GPa 0. Ref suggests that the jump in T c is also related to the /3-Po to the bcc phase 
transition and speculates that the discrepancy between the transition pressure observed 
in the room temperature X-ray experiments and the one found in the low temperature T c 
experiments is due to finite temperature effects. After the jump, T c decreases rapidly, falling 
below 4.5 K at 45 GPa |J. Enhancement in T c near structural phase transitions has been 
observed for different systems [|j], but to our knowledge no other materials exhibit such 
dramatic changes of T c with pressure. In this letter we present an ab initio calculation of 
phonon dispersion and electron-phonon interaction for the /3-Po and the bcc phases of Te. 
Our calculation reproduces and explains the experimentally observed changes in T c with 
pressure, and confirms the connection between the jump in T c and the structural phase 
transition. 

We use density functional theory (DFT) in the local density approximation (LDA). We 
compute phonon frequencies and electron-phonon coupling constants using linear-response. 
This approach has been successfully applied to the computation of phonon dispersions both 
in insulators || and in metals |||7]], and to the computation of the electron-phonon cou- 
pling ||. We employ a norm conserving pseudopotential from Ref. 0. The wavefunctions 
are expanded in a plane wave basis set with a cutoff energy of 20 Ry. Brillouin zone (BZ) in- 
tegrations are performed using Hermite-Gaussian smearing of order one 0j with a linewidth 
of 0.03 Ry and 112 (408) special k-points in the irreducible wedge of the BZ for the bcc 
(/3-Po) phase. 

From the calculated total energies we estimate the transition pressure (P c ) between the 
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/3-Po and the bcc phases to be 28 GPa, in agreement with previous LDA calculations [JTO . 
Since near the phase transition the two structures are very close, the volume-energy curves 
for these two structures are almost parallel |10] and a small change in the total energy 
for one of the phases results in a large variation of P c . Thus the error in P c due to LDA 
could be large. Indeed P c varies by 5 GPa if the total energy of one of the two phases is 
shifted by 1 mRy (a typical error in LDA). The sensitivity of P c to small energy changes 
also suggests that the difference between the observed low temperature (32-35 GPa) and 
room temperature (27 GPa) values of P c can be explained by the phonon contributions to 
the free energy. We will explore this possibility later. 

Using the approach described in Ref. 0, we compute the phonon dispersion along the 
high symmetry lines for the bcc and /3-Po phases at different volumes. For the bcc phase 
we consider four unit cell volumes, 19A 3 , 20A 3 , 21A 3 , and 22A 3 . The results are presented 
in Fig. [I]. The theoretical pressure is also shown. As the pressure decreases we observe 
an overall decrease of the phonon frequencies. In particular, the squares of the phonon 
frequencies are found to vary linearly with pressure with no peculiar behavior at P c . The 
only nearly pressure independent mode is the TN transverse phonon shown by a dot-dashed 
line. An interesting feature of the phonon dispersion is the softening of the transverse mode 
(shown by the solid line) along the TN direction. This mode exhibits a notable phonon 
anomaly, i.e. a dip, in the middle of the line. As pressure decreases, the phonon frequency 
at the dip decreases and is imaginary at a pressure of 19 GPa. This phonon anomaly is not 
related directly to the structural phase transition. Indeed, both the /3-Po and bcc phases 
have one atom per unit cell. They can not be connected with a phonon distortion with a 
non-zero q vector, since such a phonon would increase the number of atoms per unit cell. 
The fact that the phonon frequency becomes imaginary near the phase transition pressure 
is thus coincidental. Finally, at the lowest pressure another phonon anomaly along the TH 
(111) direction is visible near the F point. 

It is interesting to relate the phonon dispersion to the structure of the dynamical matrix 
in real space, D(R), R being a lattice vector. The squares of the phonon frequencies 



at a wavevector q, ^ 2 (q), are the eigenvalues of D(q)/M, M being the ionic mass and 
D(q) = -^rD(R)[1 — cos(q ■ R)]. For simplicity we restrict our analysis to the TN 
line. If q is along TN, i.e. q = pG where G is the shortest reciprocal lattice vector, and 
R is the lattice vector of a first or of a second neighbor atom, cos(q • R) is either 1 or 
cos(27rp). Thus if D(R) were short-range, i.e. D(R) ^ only for first and second neighbors, 
v 2 {p) would just be proportional to [1 — cos(27rp)] [|ll|]. The presence of a dip along TN 
indicates that the dynamical matrix also has a significant long-range contribution. In Fig. |2| 
we decompose v 2 {p) into a short-range "normal" part v 2 (p) and a long-range "anomalous" 
u aip) : u2 (p) = v nip) + v a(p) with v\{p) = A[l — cos(27rp)]. u 2 is connected to the local 
chemistry of the crystal, whereas v 2 a is a consequence of a q dependent screening induced 
by the electronic band structure. To fix the constant A in the definition of u 2 we assume 
that at the N point, i.e. the most distant point from the anomaly on the TN line, v 2 a 
is zero. Using this decomposition, Fig. |2] shows that the anomalies of the longitudinal 
mode and the transverse mode (shown by the solid line) are of equal magnitude. Moreover 
the long-range part v 2 is pressure independent whereas u 2 , exhibits a significant pressure 
dependence. Interestingly the anomaly in the phonon dispersion is more visible at low 
pressure, not because v 2 is larger but because u 2 is smaller jnfl . 

The phonon dispersion for the /3-Po structure is presented in Fig. ||] for two unit cell 
volumes, 21.92A 3 and 23A 3 . The angles a (104.90° and 104.19° respectively) are determined 
by minimizing the total energy. Some similarities can be found with the phonon dispersion 
of the bcc phase (Fig. [I]). One important difference is that for the /3-Po phase a phonon 
anomaly is present along TF but not along TL (even though both directions are equivalent 
to TN in bcc), and that the small q anomaly present in bcc on the ril (111) line is absent 
in the /3-Po on the correspondent TT line. The second important difference is the weaker 
pressure dependence of the phonon frequencies for the /3-Po phase. In particular the phonon 
frequency at the anomaly does not vary with pressure. 

Knowledge of the phonon dispersion allows us to estimate the finite-temperature free 
energy within the quasi-harmonic approximation [f[3| , F(T, V) = E(V) + F p h(T, V), where 
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E(V) is the DFT-LDA total energy and F p h(T, V) is the phonon contribution to the free 
energy: 

d?q 



FAT,V) = J2[pt kB T 
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Here i/ qj7 is the frequency of the phonon mode with wavevector q and polarization 7 com- 
puted for the unit cell volume V, and Qbz is the volume of the BZ. We use the phonon 
dispersion computed along the high symmetry lines to estimate the BZ integral fl4| . We 



then evaluate the transition pressure P c using F(T, V) instead of E(V). At T=0°K P c is 
27 GPa. At room temperature (T=293°K), P c is reduced to 23 GPa. Whereas a precise 
determination of P c is limited by the accuracy of LDA, the change in P c , AP C = 4 GPa, 
induced by the temperature is reliable. Indeed AP C is not modified when the total energy of 
one of the two phases is shifted by 1 mRy The computed value of AP C is in good agreement 
with the experimental findings. 

To study the superconducting properties we consider the electron-phonon coupling con- 



stants A q>7 [15] 
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Here e qi7 is the polarization vector, M is the atomic mass, n and m are the band indices, 
«k, n is the periodic part of the wavefunction with the eigenvalue E kjTl , N(E F ) is the density 
of states per unit cell and per spin at the Fermi energy E F , and VrV^ is the periodic part 
of the gradient of the self-consistent potential with respect to the atomic displacements with 
the wavevector q. We compute A qi7 along the high symmetry lines using Eq. (|2]), where z/ qj7 , 
e qi7 , and VrV^ are obtained from the linear- response phonon calculations. We replace the 
5-functions with Gaussians with a standard deviations of 0.021 Ry (0.035 Ry) for the bcc 
(/9-Po) phase. The BZ integration is performed using 5200 (2992) special k-points in the 
irreducible BZ of the bcc (/?-Po) structure. 

In Table | we show the computed average phonon frequency ((z^)) and the integral over 



the BZ (A) of the computed electron-phonon coupling constants |14[ . In the bcc phase A 



increases with decreasing pressure. This increase is determined completely by the variation 
of the phonon frequencies. Indeed, we find that only the prefactor l/f q 7 in Eq. ^| changes 
significantly with pressure. In the /3-Po phase A is smaller than in the bcc phase and is 
pressure independent. In particular, for the lowest phonon branches, which contribute the 
most to the value of the integrated A, both A q>7 and i/ qi7 are pressure independent. Finally 
in both phases A qi7 has peaks in q space corresponding to the phonon anomalies, in the 
bcc phase along TN and TH (111) and in the /3-Po phase along TF. In the bcc phase the 
largest contribution to A comes from the soft phonon anomaly along TO. The absence of an 
anomaly along the corresponding TT direction in the /3-Po phase is responsible for smaller 
values of A for this phase. 

Knowledge of A q 7 allows us to compute T c by numerically solving the Eliashberg equa- 
tions [|15| on a BZ mesh [14]. Our results for T c are presented in Table | and Fig. |4] for 



two typical values of the screened Coulomb interaction constant, /i*=0.12 and 0.14. The 
calculated T c is in excellent agreement with the experimental data of Ref. also shown in 
Fig. |. 

In conclusion, using DFT-LDA and linear-response, we calculate the phonon dispersion 
and the electron-phonon interaction for the bcc and the /3-Po high pressure phases of Te. 
Our results reproduce the experimental pressure dependence of the critical superconducting 
temperature, T c , and demonstrate that the observed jump in T c corresponds to the /3-Po 
to bcc structural phase transition. The pressure dependence of T c is shown to be primarily 
determined by the pressure dependence of the phonon frequencies, v. In particular: i) In /3- 
Po the lowest phonon branches, which contribute the most to the electron-phonon coupling, 
are pressure independent, so is T c . ii) The bcc phase has more soft phonon anomalies 
than /3-Po. We find an enhancement of the electron-phonon interaction at the anomalies. 
This explains the jump in T c at the phase transition, iii) In the bcc phase, the electron- 
phonon coupling varies with pressure as 1/z/ 2 . The rapid increase of the short-range part 
of the dynamical matrix with pressure results in a rapid increase of v 2 and thus in a rapid 



decrease of T c . Finally the difference of the transition pressures observed in low temperature 
and room temperature experiments is shown to arise from the phonon contribution to the 
free energy. The excellent agreement of our results with the experimental data illustrates 
the accuracy and effectiveness of DFT-LDA, linear-response approaches for study of the 
phonon-mediated superconductivity. 
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TABLES 

TABLE I. Computed average frequencies (v), electron-phonon coupling constants A, and su- 
perconducting transition temperatures T c as a function of pressure. The two values for T c corre- 
spond to two different values of fi* (0.12 and 0.14). 

P (GPa) (v) (THz) A T c (K) 



bcc 47 3.92 0.93 6.1, 5.2 

bcc 35 3.50 1.17 8.0, 7.1 

bcc 26 3.10 1.64 10.0, 9.1 

(3-Po 23 3.48 0.80 3.3, 2.7 

/3-Po 17 3.28 0.80 3.4, 2.8 
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FIGURES 

FIG. 1. The phonon dispersion for bcc. The directions of propagation and the polarizations 
are expressed in terms of the reciprocal lattice vectors. The dotted line denotes longitudinal modes. 
The solid and dot-dashed lines denote transverse modes. Along TN the polarization vector of 
the transverse mode shown by the solid line is along (111), the polarization vector of the other 
transverse mode is along (Oil). 

FIG. 2. Decomposition of v 2 into "normal" (i/^) and "anomalous" parts for two dif- 
ferent unit cell volumes of the bcc phase of Te along the TN symmetry line. The polarizations 
corresponding to different lines are given in the caption of Fig. 1. 

FIG. 3. The phonon dispersion for the /3-Po phase. The directions of propagation and the 
polarizations are expressed in terms of the reciprocal lattice vectors. Note that in the bcc lattice 
the directions (100) and (111) are equivalent to the directions (110) and (131), respectively. Along 
(110) and (111), the dotted line denotes the longitudinal modes. Along (100) and (110), the 
dot-dashed line denotes transverse modes; the polarization vector of the first of them is along 
(011), the polarization of the second depends on q. Along (110), (111), (HI), and (131), the solid 
line denotes transverse modes. The other modes do not have longitudinal or transverse character. 

FIG. 4. The critical superconducting temperature T c as a function of pressure. The exper- 
imental points Q are given by solid squares with error-bars. The computed values of T c with 
[i* = 0.12 and 0.14 are given by solid and open circles respectively. The lines are guides for the 
eye. 
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